AD-AOaa  526  UAVID  W  TAYLOR  NAVAL  SHIP  RESEARCH  ANO  DEVELOPMENT  CE— ETC  F/6  20/4 

BOUNDARY-LAYER  FLOW  ON  AN  AXISYMMETRIC  BODY  WITH  AN  INFLECTED  S— ETC(U) 
AUG  80  T  T  HUANG*  N  C  6R0VES*  G  BELT 

UNCLASSIFIED  dTNSRDC-80/064  NL 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  This  PAGE  (Wh»n  n ml*  Entered) 

~  REPORT  DOCUMENTATION  PAGE 

'  -A.EP£5I.  N.UH95.R  “  [2.  GOVT 

DTNSRDC-  86/064/  m 


2.  GOVT  ACCESSION  NO. 


mKSKDC-8fi/p6rf  fimiOifSAC 

4  TITLE J'an dSubtfff.ej  -  •••- *  -*  • 

j  BOUNDARY-LAYER  £LOW  ON  AN  ^SYMMETRIC  £ODY 
'  *  "*WITH  AN  INFLECTED  STERN  . 


7.  .AUTHORftj.  ...  _ |  .  _ _ _ -  . . . . - 

/")]  Thomas  T .  Huang {  Nancy  C.  /droves y^Garnell/Belt 


READ  INSTRUCTIONS 

_ BEFORE  COMPLETING  FORM 

'  »  “ RECIPIENT'S  Cat al og  number  " 


5.  TYPE  OF  REPORT  »  PERIOD  COVERED 


j  Final 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Bethesda,  Maryland  20084 


I  I.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 


6  PERFORMING  t)RG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUM8ERr»J 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  8  WORK  UNIT  NUMBERS 

-61152N 

Project^  ZR  GDO  M J 


(.a)  70 


U  MONITORING  AGENCY  name  ft  ADORESSflf  different  from  Controlling  Office)  IS.  SECURITY  CLASS,  (of  thie  report) 

UNCLASSIFIED 

15*7  oeCLASSlPTcATlON/ downgrading 
SCHEDULE 


ie.  distribution  statement  (of  this  Report) 

APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED 


18  supplementary  NOTES 


IS.  KEY  WORDS  (Confirm#  on  reverse  side  If  neceesmry  end  Identify  by  block  number) 

Axisymmetric  Bodies  Turbulence  Correlation  Coefficient 

Turbulent  Boundary  Layer  Turbulence  Reynolds  Stress 

(  Thick  Stern  Boundary  Layer 
^Turbulence  Length  Scale 

To  A^V^RAtT  fConttnua  on  r«v*r«»  a! dm  II  nmcmmmmry  And  Idmnllty  by  block  muitar) 

—-•Static  pressure  distributions,  mean  velocity  profiles,  distributions 
of  turbulence  intensities,  and  Reynolds  stresses,  spatial  correlation  co¬ 
efficients,  and  the  Taylor  microscale  of  turbulence  were  measured  across 
the  stern  boundary  layer  of  an  axisymmetric  body  with  an  inflected  stern. 

A  theoretical  and  numerical  iteration  scheme,  which  uses  the  boundary 
layer  and  open  wake  displacement  body,  is  found  to  model  satisfactorily  >, 


DD  ,  1473  EOITION  OF  I  NOV  «S  IS  OBSOLETE 

S^N  01 02- IF -01 4-6601 


(Continued  on  reverse  side)  j 
UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  f»hw.  Dmlm  Krlmimd) 


<  A 


7 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fHTi«n  Data  Enlarad; 


(E^ock  20  continued) 

the  viscid-inviscid  interaction  between  the  thick  stern  boundary  layer 
and  the  external  potential  flow.  The  measured  values  of  turbulence 
intensity,  eddy  viscosity,  and  mixing-length  parameters  in  the  stern 
region  are  found  to  be  much  smaller  than  those  of  a  thin  boundary 
layer.  An  approximate  similarity  characteristic  for  the  thick  axi- 
symmetric  stern  boundary  layer  is  obtained  when  the  measured  mixing- 
length  parameters,  the  measured  correlation  length  scales,  and  the 
measured  Taylor  microscales  are  normalized  by  the  square  root  of  the 
boundary-layer  cross-sectional  area  instead  of  the  local  boundary  layer 
thickness.  When  this  simple  similarity  hypothesis  for  the  mixing 
length  and  the  displacement  body  is  incorporated  into  the  McDonnell 
Douglas  Corp.,  Cebeci-Smith  differential  boundary-layer  method, 
modified  to  consider  the  displacement  body  and  wake,  the  theory  predicts 
very  well  the  measured  distributions  of  the  mean  velocity  throughout 
the  entire  stern  boundary  layer.^JThe  computation  method  developed 
earlier  was  found  to  predict  well  the  boundary  layer  on  convex  sterns 
and  is  now  found  to  apply  equally  well  to  the  flow  on  inflected  sterns. 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PACEOWlan  S«>  Bnfrtd) 


TABLE  OF  CONTENTS 


Page 

LIST  OF  FIGURES .  iv 

LIST  OF  TABLES .  v 

NOTATION .  vi 

ABSTRACT  .  1 

ADMINISTRATIVE  INFORMATION  .  1 

INTRODUCTION  .  1 

WIND  TUNNEL  AND  MODEL .  3 

INSTRUMENTATION  .  6 

COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL 

RESULTS  .  9 

MEASURED  AND  COMPUTED  PRESSURE  AND  SHEAR 

STRESS  DISTRIBUTIONS  .  10 

MEASURED  AND  COMPUTED  STATIC  PRESSURE 

DISTRIBUTION  .  12 

MEASURED  AND  COMPUTED  MEAN  VELOCITY  PROFILE  .  17 

MEASUREMENT .  17 

THEORY .  28 

COMPARISON  OF  MEAN  PROFILES .  30 

COMPARISON  OF  MEASURED  AND  COMPUTED 

INTEGRAL  PARAMETERS  .  30 

MEASURED  TURBULENCE  CHARACTERISTICS  .  34 

MEASURED  REYNOLDS  STRESSES  .  '34 

EDDY  VISCOSITY  AND  MIXING  LENGTH  .  39 

TURBULENCE  MODEL  .  41 

TURBULENCE  REYNOLDS  NUMBER  .  42 

TURBULENT  LENGTH  SCALE  IN  THICK  STERN 

BOUNDARY  LAYER  .  44 

iii 


Page 


CONCLUSIONS .  48 

ACKNOWLEDGMENT . 49 

REFERENCES .  51 


LIST  OF  FIGURES 


1  -  Three  Afterbodies  .  6 

2  -  Schematic  of  a  Two  Element  Sensor  Alined  90 

Degrees  to  Each  Other  and  45  Degrees  to 

Probe  Axis .  8 

3  -  Axisymmetric  Coordinate  Systems  .  9 

4  -  Computed  and  Measured  Stern  Pressure 

Distribution  on  Afterbody  5 .  11 

5  -  Computed  and  Measured  Shear  Stress 

Distribution  on  Afterbody  5 .  13 

6  -  Computed  and  Measured  Static  Pressure 

Distributions  Across  Stern  Boundary 

Layer  of  Afterbody  5 .  15 

7  -  Computed  and  Measured  Mean  Axial  and  Radial 

Velocity  Distribution  Across  Stern  Boundary 

Layer  of  Afterbody  5 .  31 

8  -  Comparison  of  Measured  and  Computed  Boundary 

Layer  and  Displacement  Thickness  Parameters  .  33 

9  -  Measured  Distributions  of  Reynolds  Stresses 

for  Afterbody  5 .  36 

10  -  Measured  Distributions  of  Turbulent 

Structure  Parameter  .  38 

11  -  Measured  Distributions  of  Eddy  Viscosity 

for  Afterbody  5 .  39 

12  -  Measured  Distributions  of  Mixing  Length 

for  Afterbody  5 .  40 

13  -  Similarity  Concept  for  Mixing  Length  of 

Thick  Axisymmetric  Stern  Turbulent 

Boundary  Layer  .  41 


iv 


NOTATION 


Van  Driest's  damping  factor,  A  =  26  v 


(f) 


-  1/2 


-  2 

Turbulence  structure  parameter,  a^  =  -u'v'/q 

Pressure  coefficient,  =  (P~Pq)  j  ^  pU^  =  1  -  (I^/Uq)^ 

T»  /“  o\2  T» 

Shear  stress  coefficient,  C  =  -  77-  =  7 - r- 

T  "t-P.  V" '•  /  ip 


Maximum  diameter  of  the  axisymmetric  body 

Radial  separation  distance  between  two  points  A  and  B, 

d  *  rT,  ~  r* 

r  B  A 


Total  body  length 


Bow  entrance  length 


Mixing  length  parameter  in  the  inner  region. 


£  =  0.4  Tq  £n 


fe)  l1-“p  [’  r 


Tangent  coordinate  measured  normal  to  the  body  meridian 

Value  of  n  at  which  e.  =  e_ 

1  0 

Value  of  boundary  layer  thickness  normal  'to  the 
body  surface 

Measured  local  static  pressure 
Measured  ambient  pressure 


Measured  static  pressure 


Measured  dynamic  total  pressure 


2  2  2  2 

Turbulence  parameter,  q  =  u’  +  v’  +  w’ 


Reynolds  number  based  on  model  length, 

^  X 

Axial  turbulence  Reynolds  number,  R^  =  - - — 

r  Radial  coordinate  measured  from  the  axis  of 

revolution  normal  to  the  x-axis 

Tq  Body  radius 

r  Maximum  body  radius 

max 

s  Arclength  coordinate  measured  parallel  to  the 

body  meridian 

U  Computed  potential  flow  velocity  on  the 

e  displacement  body 

Uq  Free  stream  velocity 

U -  Value  of  mean  velocity  component  in  the  axial 

J  direction  (u^)  at  the  radial  position 

u'  Turbulent  velocity  fluctuation  in  the  axial 

direction 

u  Mean  velocity  component  in  the  s  direction 

(parallel  to  the  body  meridian) 

u^  Turbulent  velocity  fluctuation  in  the  s  direction 

u'v'  Reynolds  stress 

s  n  J 

u^  Mean  velocity  component  in  the  axial  direction 

v'  Turbulent  velocity  fluctuation  in  the 

radial  direction 

v  Mean  velocity  component  in  the  n  direction 

n  (normal  to  the  body  meridian) 

v^  Velocity  fluctuation  in  the  n  direction 

vr  Mean  velocity  component  in  the  radial  direction 


*L 


vii 


Turbulent  velocity  fluctuation  in  the  azimuthal 
direction 


Axial  coordinate  measured  from  the  nose  of  the 
body  parallel  to  the  axis  of  revolution 


tan'1  (*r) 


-1 


=”  j^l+5. 5  j  ,  intermittency  factor 

Boundary  layer  thickness  where  u  /U  =  0.995, 
6  =  60.995 


Radial  boundary  layer  thickness 


*  fV  Ux 

Planar  displacement  thickness,  6=1  (1-  -rr- 

p  j  V  U6 


dr 


Axisymmetric  displacement  thickness 


Eddy  viscosity 

2  r  9us 

Eddv  viscosity  in  the  inner  region,  £  =  f — 

i  \rQJ  dn 


Eddy  viscosity  in  the  outer  region,  =  0.0168  6^ 

Nondimensional  correlation  distance 

6 


Planar  momentum  thickness,  0  = 

P 


l  (1"us)u 


dr 

6 


6r 

Axisymmetric  displacement  area,  A  =  J"  ^1-  r 


dr 


viii 


X 


Microscale  of  the  turbulence. 


u 


u 


X 


Kinematic  viscosity  of  the  fluid 
Mass  density  of  the  fluid 

Spatial  correlation  coefficient  of  two  axial  velocity 


fluctuations  at  points  A  and  B, 


pAB(dr)  =  UAUB 


Shear  stress  at  the  wall 


r  dr 


Axisyimnetric  momentum  area,  Q  = 


ABSTRACT 


Static  pressure  distributions,  mean  velocity  profiles, 
distributions  of  turbulence  intensities,  and  Reynolds 
stresses,  spatial  correlation  coefficients,  and  the  Taylor 
microscale  of  turbulence  were  measured  across  the  stern 
boundary  layer  of  an  axisymmetric  body  with  an  inflected 
stern.  A  theoretical  and  numerical  iteration  scheme, 
which  uses  the  boundary  layer  and  open  wake  displacement 
body,  is  found  to  model  satisfactorily  the  viscid-inviscid 
interaction  between  the  thick  stern  boundary  layer  and  the 
external  potential  flow.  The  measured  values  of  turbulence 
intensity,  eddy  viscosity,  and  mixing-length  parameters  in 
the  stern  region  are  found  to  be  much  smaller  than  those  of 
a  thin  boundary  layer.  An  approximate  similarity  character¬ 
istic  for  the  thick  axisymmetric  stern  boundary  layer  is 
obtained  when  the  measured  mixing-length  parameters,  the 
measured  correlation  length  scales,  and  the  measured  Taylor 
microscales  are  normalized  by  the  square  root  of  the 
boundary-layer  cross-sectional  area  instead  of  the  local 
boundary  layer  thickness.  When  this  simple  similarity 
hypothesis  for  the  mixing  length  and  the  displacement  body 
is  incorporated  into  the  McDonnell  Douglas  Corp.,  Cebeci- 
Smith  differential  boundary-layer  method,  modified  to 
consider  the  displacement  body  and  wake,  the  theory  pre¬ 
dicts  very  well  the  measured  distributions  of  the  mean 
velocity  throughout  the  entire  stern  boundary  layer.  The 
computation  method  developed  earlier  was  found  to  predict 
well  the  boundary  layer  on  convex  sterns  and  is  now  found 
to  apply  equally  well  to  the  flow  on  inflected  sterns. 


ADMINISTRATIVE  TNFORMATION 

The  work  described  in  this  report  was  funded  under  the  David  W.  Taylor 
Naval  Ship  Research  and  Development  Center's  Independent  Research  Program, 
Program  Element  61152N,  Project  Number  ZR  000  01,  and  Work  Unit  1552-103. 


INTRODUCTION 

This  report  presents  measurements  of  the  turbulent  boundary  layer 
characteristics  in  the  thick  stern  boundary  layer  of  an  axisymmetric  body 
with  an  extremely  full  inflected  stern.  These  measurements  include  the 
distributions  of  pressure  coefficient  and  shear  stress  on  the  body,  static 
pressure  distribution  across  the  stern  boundary  layer,  mean  velocity  pro¬ 
files,  turbulence  intensities,  Reynolds  stresses,  eddy  viscosity,  mixing 
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length,  and  spatial  correction  function.  This  investigation  is  an  ex¬ 
tension  of  the  earlier  research  on  convex  axisymmetr ic  thick  stern  boundary 

1  2* 

layers  reported  by  Huang,  et  al.  ’ 

The  experimental  measurements  are  compared  with  the  results  computed 

3 

by  the  method  of  Wang  and  Huang  for  computing  the  potential  flow¬ 
boundary  layer  interaction  on  axisymmetric  bodies.  The  initial  step  in 
this  procedure  begins  with  the  computation  of  the  potential-flow  pressure 

4 

distribution  on  the  body  using  the  computer  code  of  Hess  and  Smith.  The 
McDonnell  Douglas  Corp.,  Cebeci-Smith**  method"*  is  then  used  to  calculate 
the  flow  over  the  body  and  the  integral  relations  of  Granville**  are  used  to 
calculate  the  flow  in  the  wake.  The  local  displacement  thickness  computed 
in  these  boundary  layer  calculations  is  added  to  the  original  body  to 
obtain  a  modified  body  and  wake  geometry  as  suggested  by  Preston^  and 

g 

Lighthill.  The  procedure  is  repeated  until  the  pressure  distributions  on 
the  body  from  two  successive  approximations  agree  to  within  a  given  error 
criterion.  This  iterative  scheme  is  referred  to  as  the  displacement  body 
concept . 

Several  modifications  have  been  made  to  the  Douglas  C-S  differential 
5  3 

boundary  layer  method  by  Wang  and  Huang.  These  changes  were  incorporated 

to  improve  the  comparison  of  the  thick  boundary  layer  parameters  with  the 

1  9 

experimental  data  of  Huang  et  al.  and  Patel  and  Lee.  Transverse 
curvature  effects  are  included  in  the  solution  of  the  standard  thin 
boundary  layer  equations  in  the  Wang  and  Huang  method.  Effects  due  to 
longitudinal  curvature  are  neglected.  The  pressure  variation  across  the 
boundary  layer  is  obtained  from  potential  flow  calculations  for  the  final 
displacement  body.  The  mixing  length  and  eddy  viscosity  in  the  thick 
boundary  layer  region  are  obtained  by  modifying  the  thin  boundary  layer 
values  by  factors  which  are  proportional  to  the  annular  area  between  the 
body  surface  and  the  edge  of  the  boundary  layer. 

In  the  stern  and  near-wake  region  (0.95  x/L  _<  1.05),  where  x  is  the 
axial  distance  from  the  nose  and  L  is  the  total  body  length,  a  fifth  degree 

*A  complete  listing  of  references  is  given  on  page  51. 

**The  McDonnell  Douglas  Corp.,  Cebeci-Smith  method  will  hereafter  be 
designated  as  Douglas  C-S. 
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polynomial  was  used  to  determine  the  boundary  layer  flow.  The  coefficients 
for  this  polynomial  were  determined  by  the  conditions  that  the  thickness, 
slope,  and  curvature  be  equal  to  those  calculated  by  the  modified  boundary- 
layer  method  at  x/L  =  0.95  and  by  the  Granville  integral  wake  relations^  at 
x/L  =  1.05.  This  method  was  found  by  Huang  et  al.1  to  give  excellent 
agreement  with  experimental  values  of  pressure,  shear  stress,  and  velocity 
profiles  over  the  forward  90  percent  of  the  bodies  investigated.  As  the 
boundary-layer  thicknesses  became  larger,  particularly  when  these  thick¬ 
nesses  became  greater  than  the  radii  of  the  bodies,  the  measured  values  of 
shear  stress  and  velocity  became  smaller  than  those  predicted  by  the 
theory. 

In  the  following  discussion,  the  experimental  techniques  and 
geometries  of  the  model  are  given  in  detail.  The  experimental  and  theo¬ 
retical  results  are  compared  and  the  measured  turbulence  characteristics 
are  presented.  The  method  of  obtaining  the  eddy  viscosity  and  mixing 
length  is  discussed  and  the  experimentally-determined  distributions  are 
given.  The  concept  of  similarity  length-scale  in  the  thick  stern  boundary 
layer  is  examined  experimentally.  The  square  root  of  the  annular  area  be¬ 
tween  the  body  surface  and  the  edge  of  the  boundary  layer  is  found  to  be 
the  appropriate  length  scale  for  the  axisymmetric  thick  stern  boundary 
layer . 


WIND  TUNNEL  AND  MODEL 

The  experimental  investigation  was  conducted  in  the  DTNSRDC  anechoic 
wind  tunnel  facility.  The  wind  tunnel  has  a  closed  jet  test  section  that 
is  8  ft  (2.4  m)  square  and  13.75  ft  (4.19  m)  long.  The  corners  have 
fillets  which  are  carried  through  the  contraction.  The  test  section  is 
followed  by  an  acoustica±ly-lined  large  chamber  23.5  ft  (7.16  m)  long.  It 
was  found  previously  by  Huang  et  al.,^  that  the  ambient  free-stream 

turbulence  levels,  {./ u'^/Uq)  x  100,  are  0.075,  0.090,  0.100,  and  from 
0.12  to  0.15  for  free-stream  velocities  Uq  of  24.4,  30.5,  38.1,  and  45.7 
m/s,  respectively.  Integration  of  the  measured  noise  spectrum  levels  in 
the  test  section  from  10  to  10,000  Hz  indicated  that  the  typical  background 
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acoustic  noise  levels  at  30.5  m/s  were  about  93  dB  re  0.0002  dyn/cm 
(0.0002  Pa).  These  levels  of  ambient  turbulence  and  acoustic  noise  were 
considered  low  enough  so  as  not  to  unfavorably  affect  the  measurements  of 
boundary-layer  characteristics.  The  maximum  air  speed  which  can  be 
achieved  is  200  ft/sec  (61  m/s);  in  the  present  experiments  the  wind  tunnel 
velocity  was  held  constant  at  150  ft/sec  (45.72  m/s). 

An  axisymmetric  inflected  (convex  and  concave)  afterbody  with  a  bow- 
entrance  length-to-diameter  ratio  (L_/D)  of  2.05  was  used  for  the  present 

Ei 

experimental  investigation.  The  detailed  offsets  are  given  in  Table  1. 

Figure  1  shows  a  comparison  of  the  present  afterbody  (Afterbody  5)  and 

1  2 

Afterbodies  1  and  2  previously  investigated  by  Huang  et  al.  ’  The  contour 
of  Afterbody  5  changes  from  convex  to  concave  at  x/L  =  0.913,  whereas 
Afterbodies  1  and  2  are  convex  up  to  x/L  =  0.96.  Afterbody  5  is  connected 
to  a  parallel  middle  body  which  is  6.0  ft  (1.83  m)  long  and  an  existing 
streamlined  forebody  with  a  bow-entrance  length-to-diameter  ratio  (L^/D)  of 
1.82.  The  total  model  length  is  9.55  ft  (2.91  m)  with  a  maximum  forebody 
diameter  of  0.917  ft  (0.28  m) . 

The  model  was  supported  by  two  streamlined  struts  separated  by  one- 
third  of  the  model  length.  The  upstream  strut  had  a  15-cm  chord  and  the 
downstream  strut  a  3-cm  chord.  The  disturbances  generated  by  the  support¬ 
ing  struts  were  within  the  region  below  the  horizontal  centerplane,  there¬ 
fore,  all  of  the  experimental  data  were  taken  above  the  model  on  the 
vertical  centerplane  along  the  upper  meridian  where  there  was  little  effect 
from  the  supporting  struts.  One-half  of  the  model  length  protruded  beyond 
the  closed-jet  working  section  into  the  open-jet  section.  The  ambient 
static  pressure  coefficients  across  and  along  the  entire  open-jet  chamber 
(7.2  m  x  7.2  m  x  6,4  m)  were  found  to  vary  less  than  0.3  percent  of  the 
dynamic  pressure.  Tunnel  blockage  and  longitudinal  pressure  gradient 
effects  along  the  tunnel  length  were  almost  completely  removed  by  testing 
the  afterbody  in  the  open-jet  section. 

The  location  of  the  boundary-layer  transition  from  laminar  to  turbu¬ 
lent  flow  was  artificially  induced  by  a  0.024  in.  (0.61  mm)  diameter  trip 
wire  located  at  x/L  =  0.05.  Huang  et  al.^  found  that  the  trip  wire 


TABLE  1  -  MODEL  OFFSETS 


X/L  r0/L  r0/  rmax 

.0000  0.0000 

.0105  .2t93 

.0150  .3118 

.0184  .3835 

.0213  .4441 

.0239  .4975 

.0262  .5454 

.0283  .5891 

.0302  .6291 

.0320  .6659 

.0336  .7000 

.0351  .7315 

.0365  .7607 

.0378  .7877 

.0390  .8126 
.0401  .8356 

.0411  .8567 

.0420  .8760 

.0429  .8936 

.0437  .9097 

.0443  .9241 

.0450  .9371 

.0455  .9486 

.0460  .9587 

.0464  .9675 

.0468  -9251 

.0471  .9816 

.0474  .9869 

.0476  .9911 

.0477  .9945 

.0478  .9969 

.0479  .9986 

.0480  .9996 

.0480  I ■ 0000 
.0480  1 • 0000 

.0480  1 • 0000 

.0480  1 .0000 

.0480  1 .0000 

. 0480  1 . 0000 

.0480  1 . 0000 

.0480  1 .0000 

.0480  1 .0000 

. 0480  1 . 0000 

.0480  I .0000 
.0480  I .0000 
.0480  1 .0000 
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Figure  1  -  Three  Afterbodies 

effectively  tripped  the  flow  0.39  in.  (1  cm)  downstream  from  the  wire  and, 
because  of  the  finite  parasitic  drag  of  the  wire,  the  boundary-layer  can  be 
considered  to  become  turbulent  at  a  virtual  origin  5.4  in.  (13.7  cm)  up¬ 
stream  from  the  trip  wire  at  a  Reynolds  number  of  9.3  *  10^.  The  virtual 
origin  for  the  turbulent  flow  is  defined  such  that  the  sum  of  the  laminar 
frictional  drag  from  the  nose  to  the  trip  wire,  the  parasitic  drag  of  the 
trip  wire,  and  the  turbulent  frictional  drag  aft  of  the  trip  wire  is  equal 
to  the  sum  of  the  laminar  frictional  drag  from  the  nose  to  the  virtual 
origin  and  the  turbulent  frictional  drag  from  the  virtual  origin  to  the 
after  end  of  the  model  (similar  to  Reference  10) .  The  virtual  origin 
location  was  taken  as  the  transition  location  at  which  the  boundary-layer 
changed  from  laminar  to  turbulent  in  the  boundary-layer  calculations. 

INSTRUMENTATION 

A  Preston  tube  with  a  0.072  in.  (1.83  mm)  inside  diameter  was  attached 
and  alined  with  the  flow  on  the  upper  meridian  of  the  stern  to  measure  the 
shear  stress  distribution.  The  Preston  tube  was  calibrated  in  a  1  in. 

(2.54  cm)  diameter  water-pipe  flow  facility  described  by  Huang  and  von 
Kerczek.11  A  series  of  0.031  in.  (0.8  mm)  diameter  pressure  taps  were 
embedded  on  the  upper  meridian  of  the  stern  at  the  Preston  tube  locations. 
These  pressure  taps  were  connected  to  a  multiple  pressure  scanivalve  system 
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that  takes  one  integral  pressure  transducer  with  its  zeroing  circuit  and 

measures  a  single  pressure  in  sequence  along  the  stern  upper  meridian. 

The  pressure  transducer  used  was  designed  for  measuring  low  pressure  up  to 
3 

0.2  psi  (1.379  x  10  Pa).  The  zero-drift  linearity,  scanivalve,  hystere¬ 
sis,  and  pressure  transducer  with  zeroing  circuit  were  carefully  checked 
and  the  overall  accuracy  was  found  to  be  within  0.5  percent  of  the  dynamic 
pressure. 

A  Prandtl-type  static  pressure  probe  of  0.123  in.  (3.125  mm)  diameter 
with  four  equally  spaced  holes  located  three  probe  diameters  after  the 
probe  nose  was  used  to  measure  the  static  pressure  across  the  boundary 
layer.  By  yawing  the  probe  in  the  free  stream,  it  was  found  that  the 
measured  static  pressure  was  insensitive  to  the  probe  angle  up  to  a  5- 
degree  yaw.  In  all  cases,  the  static  pressure  probe  remained  in  alinement 
by  less  than  5-degree  yaw. 

The  mean  axial  and  radial  velocities  and  the  turbulence  intensities 
for  the  Reynolds  stress  calculations  were  measured  by  a  TSI,  Inc.  Model 
1241-20  "X"  type  hot-film  probe.  The  probe  elements  are  0.002  in.  (0.05 
mm)  in  diameter  with  a  sensing  length  of  0.04  in.  (1.0  mm).  The  spacing 
between  the  two  cross  elements  is  0.04  in.  (1.0  mm).  A  typical  schematic 
of  the  hot-film  probe  used  is  shown  in  Figure  2.  A  two-channel  hot-wire 
and  hot  film  anemometer  with  linear izers  was  used  to  monitor  the  response 
of  the  hot-film  probe.  A  temperature  compensating  sensor  (probe)  was  used 
with  each  hot-film  element  to  regulate  the  operating  temperature  of  the 
sensor  with  the  changes  in  air  temperature.  The  "X"  hot  film  and  its 
temperature-compensated  sensor  were  calibrated  together  through  the  ex¬ 
pected  air  temperature  range  and  supplied  with  their  individual  lineariza¬ 
tion  polynomial  coefficients  at  the  factory.  A  single  element  TSI,  Inc. 
Model  1212-TI.5  hot-wire  probe  also  was  used  to  measure  the  velocity 
distribution  across  the  stern  boundary  layer.  Two  single-element  TSI 
(Model  1212)  hot-wire  probes  with  a  90-deg  bend  separated  by  a  small  radial 
distance  also  were  used  to  measure  spatial  correlation  functions  in  the 
thick  boundary  layer. 
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PERPENDICULAR  TO  MEAN 
FLOW  =  -Jv* 

FILM  SENSOR:  0  002  IN.  DIA  (0.05  mm) 

Figure  2  -  Schematic  of  a  Two  Element  Sensor  Alined  90  Degrees 
to  Each  Other  and  45  Degrees  to  Probe  Axis 

The  frequency  response  of  the  anemometer  system,  for  reliable  measure 
ments  claimed  by  the  manufacturer,  is  0  to  200  kHz,  Calibration  of  the  "X 

hot  film  was  made  before  and  after  each  set  of  measurements.  It  was  found 

that  the  hot-film  anemometer  system  had  a  +0.5  percent  accuracy,  +0.75 

ft/sec  (+0.23  m/s)  accuracy  at  the  free  stream  velocity  of  150  ft/sec 

(45.72  m/s),  during  the  entire  experiment.  An  estimate  was  made  of  the 
cross-flow  velocity  by  yawing  the  "X"  hot-film  probe  in  the  freestream.  1 
was  found  that  the  crossflow  velocities  were  about  one  percent  of  the 
freestream  velocity. 

The  linearized  signals  were  fed  into  a  Time/Data  Model  1923-C  real¬ 
time  analyzer.  Both  channels  of  the  analog  signal  were  digitized  at  a 
rate  of  128  points  per  second  for  8  seconds.  These  data  were  immediately 
analyzed  by  a  computer  to  obtain  the  individual  components  of  mean 
velocity,  turbulence  fluctuation,  and  Reynolds  stress  on  a  real  time  basis 
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A  traversing  system  with  a  streamlined  strut  was  mounted  on  a  guide 
plate  that  permitted  the  traverse  to  be  locked  in  various  stationary 
positions  parallel  to  the  longitudinal  model  axis. 

COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  RESULTS 
Two  coordinate  systems  for  an  axisymmetric  body  are  given  in  Figure  3. 
The  coordinate  system  used  to  experimentally  measure  the  boundary-layer 


Figure  3  -  Axisymmetric  Coordinate  Systems 


flow  and  to  report  the  results  is  the  x-r  coordinate  system.  The  axial 

coordinate  x  is  measured  from  the  nose  of  the  body  parallel  to  the  axis  of 

revolution.  The  radial  component  r  is  measured  from  the  axis  of  revolution 

and  is  normal  to  the  x-axis.  The  curvilinear  s-n  coordinate  system  is  used 

in  the  Douglas  C-S  differential  boundary  layer  method. The  arclength 

coordinate  s  is  measured  parallel  to  the  body  meridian  and  the  tangent 

coordinate  n  is  measured  normal  to  the  body  meridian. 

The  Douglas  C-S  boundary  layer  method  is  used  in  conjunction  with  the 

displacement  body  concept  to  predict  the  flow  over  the  body.  Wang  and 
3 

Huang  have  made  several  modifications  to  the  Douglas  C-S  code  in  an 
attempt  to  better  model  a  thick  stern  boundary  layer.  A  comparison  of 
their  original  theory  with  experimental  results  is  given  by  Huang  et  al.^ 
That  comparison  revealed  several  inadequacies  in  the  theory.  Wang  and 

3 

Huang  subsequently  modified  their  original  method  based  on  these  experi¬ 
mental  comparisons.  The  new  modifications  are  evaluated  in  the  present 
investigation.  The  mixing  length  in  the  thick  stern  boundary  layer  region 
is  modeled  by  a  simple  algebraic  formula  proposed  by  Huang  et  al.^  This 
formulation  includes  only  transverse  curvature  effects.  The  eddy  viscosity 


in  the  outer  region  was  also  modified  by  an  algebraic  formula  with  the  in¬ 
tention  of  improving  the  accuracy  of  the  velocity  profiles  near  the  body 
in  the  thick  boundary  layer  region.  The  final  modification  to  the  original 
method  uses  the  predicted  potential  flow  velocities  outside  of  the  dis¬ 
placement  body  to  modify  the  velocities  in  the  direction  parallel  to  the 
body  meridian  tangent. 


MEASURED  AND  COMPUTED  PRESSURE  AND  SHEAR 
STRESS  DISTRIBUTIONS 


Pressure  taps  were  used  to  measure  the  steady  pressure  at  15  selected 
locations  along  the  stern  surface  of  the  axisymmetric  body.  The  pressure 
coefficient,  C  ,  is  computed  from  the  measured  pressures  by  the  relation¬ 
ship 


C 

P 


P-po 

Pt-ps 


P“Pr 


2  P 


(1) 


where  p  =  measured  local  static  pressure 
pQ  =  measured  ambient  pressure 
Pt  =  measured  dynamic  total  pressure 
pg  =  measured  static  pressure 
p  =  mass  density  of  the  fluid 

The  analytically-predicted  pressure  coefficient  was  computed  on  the  dis¬ 
placement  body  using  the  Douglas-Neumann  potential  flow  method  of  Hess  and 
4 

Smith  and  is  given  by 


where  U^  is  the  computed  potential  flow  velocity  on  the  displacement  body 
and  UQ  is  the  freestream  velocity,  150  ft/sec  (45.72  m/s). 
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A  comparison  of  the  measured  and  theoretical  values  of  the  pressure 
coefficient  on  the  surface  of  the  axisymmetric  body  under  consideration  is 
given  in  Figure  4.  The  theory  underpredicts  the  magnitude  of  the  minimum 


Figure  4  -  Computed  and  Measured  Stern  Pressure  Distribution  on 

Afterbody  5 


trough.  The  three  measured  values  of  C  for  •=-  >  0.975  fluctuate  somewhat 

p  L 

about  the  theoretical  curve  for  C^.  These  discrepancies  are  small,  and 
overall  agreement  between  theoretical  and  measured  values  of  the  pressure 
coefficient  is  considered  excellent. 

Readings  from  a  Preston  tube,  which  was  taped  to  the  stern  surface  at 
the  pressure  tap  locations,  were  used  in  conjunction  with  the  steady 
pressure  readings  to  obtain  the  shear  stress  distribution  at  the  body  sur¬ 
face.  The  calibration  curve,  presented  by  Huang  and  von  Kerczek,^  for  a 


Preston  tube  in  a  pipe  was  used  to  compute  the  shear  stress  at  the  wall, 

T  .  The  shear  stress  coefficient  C  is  given  by 
W  T 


C 

I 


Tw  /U0 
Pt-ps\Ue 


(3) 


for  the  measured  shear  stress  and  by 


C 

T 


(4) 


for  the  analytical  shear  stress. 

Figure  5  shows  a  comparison  of  the  measured  and  theoretical  values  of 
the  shear  stress  coefficient.  Agreement  is  good  for  values  of  —  <  0.9. 

The  trough,  which  occurs  for  —  >  0.9,  is  somewhat  overpredicted  by  the 
theory.  This  discrepancy  in  the  distribution  of  the  shear  stress  co¬ 
efficients  indicates  that  the  present  analytical  model  may  be  inadequate 
for  the  precise  prediction  of  the  shear  stress  on  the  last  ten  percent  of 
the  body.  The  measured  values  of  C  and  C  are  tabulated  in  Table  2. 

J  p  T 

MEASURED  AND  COMPUTED  STATIC  PRESSURE 
DISTRIBUTION 

The  measured  and  computed  static  pressure  coefficients  for  Afterbody  5 
are  compared  in  Figure  6  at  various  locations  across  the  stern  boundary 
layer.  The  off-body  option  of  the  Douglas  potential-flow  computer  code  was 
used  to  compute  the  static  pressure  distributions  for  the  displacement  body 
(solid  lines)  and  the  actual  body  (broken  lines) .  As  can  be  seen  in 
Figure  6,  except  at  x/L  =  0.704,  0.987,  and  1.045,  the  measured  static 
pressure  distributions  agree  better  with  the  theoretical  pressure  distri¬ 
butions  computed  from  the  displacement-body  model  than  those  from  the 
original  body.  At  x/L  =  0.704,  0.987,  and  1.045,  the  two  calculation 
methods  agree  to  within  one  percent.  The  discrepancy  between  the  measured 
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L 

Figure  5  -  Computed  and  Measured  Shear  Stress  Distribution 

on  Afterbody  5 


13 


TABLE  2  -  MEASURED  PRESSURE  AND  SHEAR  STRESS  COEFFICIENTS 


x/L 

rn/r 

0  max 

C 

P 

C 

T 

0.7036 

1.0000 

-0.024 

0.00268 

0.7507 

1.0000 

-0.048 

0.00272 

0.7768 

1.0000 

-0.075 

0.00278 

0.8000 

1.0000 

-0.140 

0.00294 

0.8135 

0.9975 

-0.200 

0.00300 

0.8308 

0.9756 

-0.238 

0.00295 

0.8465 

0.9296 

-0.220 

0.00274 

0.8727 

0.7983 

-0.070 

0.00188 

0.8884 

0.6939 

+0.030 

0.00174 

0.9093 

0. 5378 

+0.133 

0. 00069 

0.9303 

0.3784 

+0.182 

0.00043 

0.9512 

0.2345 

+0.201 

0.00034 

0.9774 

0.1199 

+0.186 

0.00058 

0.9874 

0.1108 

+0.146 

0.00101 

1.0000 

0.0000 

+0.173 

— 
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Figure  6  -  Computed  and  Measured  Static  Pressure  Distributions  Across  Stern 

Boundary  Layer  of  Afterbody  5 


Figure  6a  -  Nondimensional  Axial  Lengths,  x/L  =  0.704, 


0.831,  0.873 
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Figure  6  (Continued) 


Cp 

Figure  6b  -  Nondimensional  Axial  Lengths,  x/L  =  0.951,  0.909 
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Figure  6c  -  Nondimensional  Axial  Lengths,  x/L  =  0.987,  1.045 
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and  computed  (by  displacement-body  model)  values  of  is,  in  general,  less 

than  0.01  which  is  about  the  accuracy  of  the  measurement.  The  measured 

values  of  C  are  tabulated  in  Table  3. 

P 

As  will  be  seen  later,  the  displacement  body  of  Afterbody  5  is  convex 
from  the  parallel  middle  body  up  to  x/L  =  0.898  and  becomes  concave  down¬ 
stream  from  x/L  >  0.90,  which  is  rather  close  to  the  inflection  point  of 
the  actual  body  (x/L  =  0.913).  The  differences  in  the  theoretical  values 
of  off-body  pressure  distributions  computed  from  the  displacement  body  and 
the  actual  body  are  very  noticeable  at  x/L  =  0.951  where  the  measured 
pressure  distributions  agree  rather  well  with  the  distributions  computed 
from  the  displacement-body  model.  The  close  agreement  between  the  computed 
displacement-body  and  measured  static  pressure  distributions  for  Afterbody 
5,  as  well  as  Afterbodies  1  and  2  (Reference  1),  support  the  displacement- 
body  concept  for  computing  the  potential  flow  outside  of  the  displacement 
surface.  The  pressure  coefficients  inside  the  displacement  body  may  be 

more  accurately  calculated  by  a  streamline  curvature  method  proposed  by 
12 

Dyne.  However,  in  the  present  calculation  method  the  value  of  C  is 
assumed  to  be  constant  inside  the  displacement  body.  The  maximum  error  of 
the  value  of  C  from  this  approximation  is  less  than  2  percent  inside  the 
displacement  body  (Figure  6). 

MEASURED  AND  COMPUTED  MEAN  VELOCITY  PROFILE 

MEASUREMENT 

The  axial  and  radial  mean  velocity  components  were  measured  by  the 

"X"  hot-film  sensor.  At  a  given  axial  location,  measurements  were  taken 

at  numerous  radial  positions  to  obtain  the  velocity  profile.  The  hot-film 

probe  was  mounted  with  the  sensing  elements  in  a  vertical  position.  At 

each  radial  position,  1024  data  values  were  collected  in  8  seconds.  The 

linearized  data  were  input  to  a  computer  code  which  calculated  and  printed 

the  values  of  the  mean  velocity  components  u^  in  the  axial  direction  and 

v  in  the  radial  direction  (given  in  Table  4). 
r 

A  single  element  hot-wire  sensor  also  was  used  to  compare  with  the 
results  of  the  "X"  sensor.  The  values  of  the  axial  mean  velocity  component 
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TABLE  3  -  MEASURED  STATIC  PRESSURE  COEFFICIENTS  ACROSS  STERN 
BOUNDARY  LAYER  AND  NEAR  WAKE 


x/L  =  0.7036  x/L  =  0.8308 


r/r 

max 

(r-r_)/r 

0  max 

C 

P 

r/r 

max  ! 

(r-rn)/r 

0  max 

C 

P 

1.0000 

0.0000 

-0.0240 

0.9756 

0.0000 

-0.2380 

1.0873 

0.0873 

-0.0242 

0.9938 

0.0182 

-0.2372 

1.1576 

0.1576 

-0.0251 

1.0711 

0.0955 

-0.2072 

1 .2349 

0.2349 

-0.0222 

1.1484 

0.1728 

-0.1821 

1.2981 

0.2981 

-0.0242 

1.2187 

0.2431 

-0.1599 

1.3614 

0.3614 

-0.0232 

1.2890 

0.3134 

-0.1412 

1.4317 

0.4317 

-0.0242 

1.3663 

0.3907 

-0.1240 

1.5090 

0.5090 

-0.0241 

1.4366 

0.4610 

-0.1132 

1.5792 

0.5792 

-0.0241 

1.4998 

0.5242 

-0.1019 

1.6495 

0.6495 

-0.0227 

1.5631 

0.5875 

-0.0901 

1.7198 

0.7198 

-0.0232 

1.6404 

0.6648 

-0.0827 

1.8252 

0.8252 

-0.0232 

1.7458 

i 

0.7702 

-0.0723 

1.9377 

0.9377 

-0.0227 

1.8442 

0.8686 

-0.0625 

2.0642 

1.0642 

-0.0207 

1.9567 

0.9810 

-0.0551 

2.2118 

1.2118 

-0.0217 

2.0621 

1.0865 

-0.0482 

2.3523 

1.3523 

-0.0197 

2.2027 

1.2270 

-0.0433 

2.4929 

1.4929 

-0.0207 

2.3432 

1.3676 

-0.0354 

2.6335 

1.6335 

-0.0182 

2.4978 

1.5222 

-0.0300 

2.7740 

1.7740 

-0.0172 

2.6595 

1.6839 

-0.0251 

2.9497 

1.9497 

-0.0168 

2.8352 

1.8596 

-0.0222 

3.1254 

2.1254 

-0.0138 

3.1233 

2.1477 

-0.0128 

3.3082 

2.3082 

-0.0153 

3.3834 

2.4078 

-0.0103 

3.4768 

2.4768 

-0.0133 

3.6083 

2.6327 

-0.0049 

3.7228 

2.7228 

-0.0123 

TABLE  3  (Continued) 


x/L  =  0.8727 

x/L  =  0.9093 

r/r 

(r-rn)/r 

C 

r/r 

(r-r _) /r 

c 

max 

0  max 

P 

max 

0  max 

P 

0.7983 

0.0000 

-0.0700 

0.5378 

0.0000 

+0.1330 

0.8220 

0.0236 

-0.0709 

0.5687 

0.0309 

+0.1318 

0.8853 

0.0869 

-0.0680 

0.6390 

0.1012 

+0.1190 

0.9555 

0.1572 

-0.0611 

0.7093 

0.1715 

+0.1067 

1.0188 

0.2204 

-0.0527 

0.7866 

0.2488 

+0.0984 

1.0961 

0.2977 

-0.0480 

0.8428 

0.3050 

+0.0920 

1.1594 

0.3610 

-0.0434 

0.8780 

0.3402 

+0.0984 

1.2296 

0.4313 

-0.0409 

0.9272 

0.3893 

+0.0831 

1.2999 

0.5016 

-0.0365 

0.9834 

0.4456 

+0.0772 

1.3772 

0.5789 

-0.0350 

1.0466 

0.5088 

+0.0718 

1.4475 

0.6491 

-0.0301 

1.0958 

0.5580 

+0.0669 

1.5178 

0.7194 

-0.0287 

1.1731 

0.6353 

+0.0610 

1.6162 

0.8178 

-0.0252 

1.2294 

0.6916 

+0.0566 

1.7216 

0.9232 

-0.0192 

1.3840 

0.8462 

+0.0442 

1.8622 

1.0638 

-0.0183 

1.5175 

0.9797 

+0.0344 

2.0027 

1.2044 

-0.0148 

1.7003 

1.1624 

+0.0255 

2.1433 

1.3449 

-0.0104 

1.8619 

1.3241 

+0.0206 

2.4525 

1.6542 

-0.0074 

2.0446 

1.5068 

+0.0152 

2.8110 

2.0126 

-0.0044 

2.2203 

1.6825 

+0.0132 

3.1624 

2.3640 

+0.0010 

2.3960 

1.8582 

+0.0097 

3.4505 

2.6522 

+0.0049 

2.7475 

2.2096 

+0.0072 

2.9583 

2.4205 

+0.0062 

3.2043 

2.6665 

+0.0067 
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TABLE  3  (Continued) 


x/L  =  0.9512 

x/L  =  0.9874 

r/r 

(r-r„)/r 

C 

r/r 

(r-rn)/r 

C 

max 

0  max 

P 

max 

0  max  i 

j 

P 

0.2345 

0.0000 

+0.2010 

0.1108 

0.0000 

+0.1460 

0.2655 

0.0309 

+0.1845 

0.1245 

0.0136 

+0.1488 

0.2795 

0.0450 

+0.1869 

0.1527 

0.0418 

+0.1497 

0.3287 

0.0942 

+0.1802 

0.1737 

0.0628 

+0.1507 

0.3920 

0.1574 

+0.1751 

0.2300 

0.1191 

+0.1536 

0.4693 

0.2347 

+0.1642 

0.2721 

0.1612 

+0.1517 

0.5325 

0.2980 

+0.1549 

0.3494 

0.2385 

+0.1497 

0.6028 

0.3683 

+0.1430 

0.4197 

0.3088 

+0.1448 

0.6801 

0.4456 

+0.1351 

0.4900 

0.3791 

+0.1388 

0.7434 

0.5088 

+0.1243 

0.5603 

0.4494 

+0.1304 

0.8136 

0.5791 

+0.1144 

0.6306 

0.5197 

+0.1225 

0.9542 

0.7197 

+0.1011 

0.7009 

0.5899 

+0.1171 

1.0948 

0.8602 

+0.0848 

0.7711 

0.6602 

+0.1092 

1.2353 

1 . 0008 

+0.0719 

0.8414 

0.7305 

+0.0978 

1.4462 

1.2116 

+0.0601 

0.9539 

0.8430 

+0.0904 

1.6641 

1.4295 

+0.0503 

1.0241 

0.9132 

+0.0865 

1.8679 

1.6333 

+0.0429 

1.2350 

1.1241 

+0.0727 

2.0857 

1.8512 

+0.0350 

1.3685 

1.2576 

+0.0678 

2.2966 

2.0620 

+0.0310 

1.5864 

1.4755 

+0.0555 

2.5004 

2.2659 

+0.0256 

1.7902 

1.6793 

+0.0481 

2.7112 

2.4767 

+0.0241 

2.0011 

1.8902 

+0.0422 

2.9291 

2.6946 

+0.0212 

2.2119 

2.1010 

+0.0387 

3.1400 

2.9054 

+0.0202 

2.4298 

2.3189 

+0.0353 

2.6336 

2.5227 

+0.0304 

2.8444 

2.7335 

+0.0264 

3.1958 

3.0849 

+0.0230 

3.4770 

3.3661 

+0.0210 
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TABLE  3  (Continued) 


x/L  = 

1.045 

r/r 

C 

max 

P 

0.0000 

+0.0404 

0.0562 

+0.0444 

0.1195 

+0.0483 

0.1687 

+0.0518 

0.2179 

+0.0537 

0.2811 

+0.0562 

0.3374 

+0.0557 

0.4076 

+0.0552 

0.4849 

+0.0562 

0.5552 

+0.0547 

0.6255 

+0.0547 

0.6888 

+0.0517 

0.7731 

+0.0497 

0.8715 

+0.0472 

0.9699 

+0.0467 

1.0823 

+0.0427 

1.1807 

+0.0427 

1.2932 

+0.0387 

1.3986 

+0.0382 

1.5040 

+0.0370 

1.6516 

+0.0323 

1.8625 

+0.0298 

2.0311 

+0.0288 

2.3474 

+0.0253 

2.6356 

+0.0209 

t 


TABLE  4  -  MEASURED  MEAN  AND  TURBULENT  VELOCITY  CHARACTERISTICS 


(Continued) 


(Continued) 


TA3LE  4  (Continued) 


TABLE  4  (Continued) 


x/L  -  0.951 


0.2346  tan  a  ■  -  0.2952 


7  u  7  V  r  w 


r-ro 

£ 

£ 

100  'u'v' 

— u  'v ' 

uo 

2 

q 

6 

r 

* 

M 

6  p 

6 

r 

/ <r0+6r)2-r 

0.0443 

0.154 

0.0347 

0.0004 

0.0054 

0.0041 

0.0584 

0.167 

0.0657 

0.0011 

0.0127 

0.0098 

0.0775 

0.191 

0.0864 

0.0015 

0.0178 

0.0137 

0.0852 

0.198 

0.1018 

0.0024 

0.0219 

0.0169 

0.0830 

0.178 

0.1277 

0.0026 

0.0245 

0.0188 

0.0924 

0.189 

0.1482 

0.0026 

0.0226 

0.0174 

0.1053 

0.194 

0.1792 

0.0027 

0.0227 

0.0174 

0.0996 

0.177 

0.2102 

0.0026 

0.0223 

0.0172 

0.1077 

0.190 

0.2515 

0.0029 

0.0235 

0.0181 

0.0973 

0.182 

0.3135 

0.0027 

0.0234 

0.0180 

0.0744 

0.161 

0.3342 

0.0024 

0.0236 

0.0182 

0.0729 

0.161 

0.3703 

0.0035 

0.0349 

0.0269 

0.0791 

0.179 

0.4063 

0.0038 

0.0366 

0.0282 

0.0664 

0.160 

0.4373 

0.0036 

0.0369 

0.0284 

0.0672 

0.170 

0.4890 

0.0041 

0.0423 

0.0326 

0.0568 

0.167 

0.5870 

0.0040 

0.0446 

0.0343 

0.0436 

0.175 

0.6954 

0.0027 

0.0346 

0.0268 

0.0130 

0.130 

0.8090 

0.0014 

0.0331 

0.0255 

0.0001 

0.003 

0.9432 

0.0000 

0.0100 

0.0078 

^  =  0.0480. 


TABLE  4  (Continued) 


0.0182 
I  0.0252 
0.0322 
0.0513 
,  0.0604 
.  0.0814 
0.1025 
0.1306 
0. 1482 
1  0.1763 

0.2290 
0.2853 
0. 3520 
0.4610 
0.5383 
0.6085 
0.6823 
0.  7456 
0.8124 
0.8897 
1.0583 
1.2305 
1 . 5890 
1.9369 
2 . 3094 
2.6397 
2.8013 


j  0. 39^ 

-0.021 

0.412 

-0.025 

0 . 4  34  j 

-0.025 

i  0.465 

-0.026 

.  0.487 

1  -0.028 

1  0.520 

-0.028 

1  0. 593 

1  -0.030 

j  0.655 

:  -0.031 

0.721 

'  -0.030  J 

0.809  , 

1  -0.031  J 

0.855 

-0.029  1 

0.899 

1  -0.028  | 

0.929 

-0.025 

0.938 

-0.025 

0.941 

-0.022 

0.946 

-0.020 

0.955  i 

-0.017 

0.959  | 

-0.013 

0.961 

i  -0.006 

0.971 

-0.001 

0.977  | 

0.001 

0.977 

0.008 

0.979 

0.010 

1 

measured  by  the  single-element  probe  are  presented  in  the  figures  along 
with  the  values  measured  by  the  two-component  probe. 


THEORY 

The  Douglas  C-S  method^  solves  the  standard  incompressible,  steady 
continuity  and  momentum  equations  for  thin  axisymmetr ic  boundary  layers. 
This  set  of  partial  differential  equations,  with  the  addition  of  the 
transverse  curvature  effects,  is 


3(ru  )  3(rv  ) 

- s_  +  - 

3s  3n 


0 


(5) 


and 


u 

s 


+  v 

n 


1  d£  +  1  3_ 
p  ds  r  3n 


(6) 


where  u 

s 


v 


n 


P 

r 

r 

X 

V 


=  mean  velocity  component  in  the  s  direction  (parallel 
to  the  body  meridian) 

=  mean  velocity  component  in  the  n  direction  (normal  to 
the  body  meridian) 

=  pressure  on  the  body 

=  rQ(s,n)  +  n  cos  a 

=  body  radius 

=  axial  distance  measured  from  the  nose 
=  kinematic  viscosity  of  the  fluid 


u'v'  =  Reynolds  stress 
s  n 

u^  =  velocity  fluctuation  in  the  s  direction 
V  =  velocity  fluctuation  in  the  n  direction 
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1 


The  wall  boundary  conditions  are  u^  (s,n)  =  (s,n)  =  0  at  n  =  0.  Effects 

due  to  longitudinal  curvature  are  neglected.  The  Reynolds  stress  model 
is  given  in  detail  in  the  section  presenting  the  turbulence  characteristics. 

The  Reynolds  stress  -u'v'  in  the  modified  Douglas  C-S  computer  code 

1  5  13  S  n 

is  modeled  ’’  by 


s  n 


£ .  ' 
l 


u 


3u 

c 

3n 


for  0  <  n  <  n 
—  —  c 


for  n  <  n  <  6 
c  —  — 


(7) 


where  £ 


^  ft) 


3u 

c 

3n 


eddy  viscosity  in  the 
inner  region 


e.  =  0.0168  U  6  Y 
0  e  p  'tr 


eddy  viscosity  in  the 
outer  region 


mixing  length  parame¬ 
ter  in  the  inner  region 


Van  Dr Test's  damping 
factor 

boundary  layer  thick- 
u 

ness  ^  =  0. 995 
e 


planar  displacement 
thickness 

computed  potential  flow 
velocity  at  the  edge 
of  the  displacement 
body 

wall  shear  stress 
value  of  n  at  which 
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The  Douglas  C-S  computer  code  is  used  to  calculate  the  velocity 
components  ug  and  parallel  to  and  normal  to  the  body  surface,  respec¬ 
tively.  Wang  and  Huang^  modified  these  velocity  components  by  using  the 
predicted  potential  flow  velocities  outside  of  the  displacement  body  and 
converted  the  normal  and  tangential  velocity  components  to  axial  and  radial 
components  for  comparison  with  the  experimental  values. 

COMPARISON  OF  MEAN  PROFILES 

Figure  7  shows  a  comparison  of  the  theoretically  predicted  and 

experimentally  measured  mean  axial  u^  and  radial  v  velocity  profiles  at 

various  axial  locations  on  the  body.  The  velocity  components  are  non- 

dimensionalized  by  the  free  stream  velocity  U^.  The  solid  curves  represent 

the  theoretical  results,  the  open  symbols  are  the  "X"  hot-film  values,  and 

the  solid  symbols  are  the  single  element  hot-wire  values.  The  accuracies 

of  the  experimental  measurements  of  u  /IL  and  v  /U_  were  estimated  to  be 
r  x  0  r  0 

about  0.5  percent  and  1.0  percent,  respectively. 

Figure  7  shows  the  good  agreement  between  the  theoretical  and  the 

"X"  hot-film  experimental  axial  velocity  profiles,  and  even  better  agreement 

3 

with  the  theory  for  the  single  element  sensor.  Wang  and  Huang  modified 
the  eddy  viscosity  model,  and  this  modeling  plays  an  important  role  in 
computations  of  axial  velocity  profiles  near  the  wall.  The  agreement  be¬ 
tween  theory  and  experiment  is  shown  in  Figure  7  to  be  very  satisfactory 
near  the  wall. 

The  agreement  between  theory  and  experiment  for  the  radial  velocity 
profiles  is  generally  less  satisfactory.  This  velocity  component  is  small 
and  more  difficult  to  measure  accurately.  Nevertheless,  the  maximum 
difference  between  the  computed  and  measured  values  of  v^/Uq  is  less  than 
0.03. 

COMPARISON  OF  MEASURED  AND  COMPUTED  INTEGRAL  PARAMETERS 

& 

The  distributions  of  (radial)  displacement  thickness  6  and  (radial) 
boundary  layer  thickness  6  are  given  in  Figure  8  along  with  the  stern 
profile  of  the  axisymmetric  body.  The  theoretical  boundary  layer  thickness 
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Figure  7a  -  Nond imensional  Axial  Lengths,  x/L  =  0.704,  0.831,  0.873 


Nondimensional  Axial  Lengths,  x/L  =  0,909, 


L 

Figure  8  -  Comparison  of  Measured  and  Computed  Boundary  Layer 
and  Displacement  Thickness  Parameters 


is  defined  as  the  normal  distance  from  the  body  where  the  nondimen sionalized 
tangential  velocity  u  /U  has  a  value  of  0.995.  In  this  definition,  u  is 

S  6  S 

the  dimensional  tangential  velocity  and  Ug  is  the  dimensional  local 
tangential  velocity  calculated  by  the  potential  flow  program  using  the 
latest  body-wake  displacement  model. 

The  measured  value  of  the  boundary  layer  thickness  6^  is  a  difficult 
value  to  obtain.  This  variable  is  defined  as  the  radial  position  where 
Reynolds  stress  uN/^  becomes  "small."  The  overall  accuracy  of  the  measured 
values  of  6^  presented  in  the  present  investigation  is  estimated  to  be 
about  10  percent. 

The  integral  parameters  of  interest  are:*- 
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where  U.  is  the  value  of  u  at  r  =  r.  +  5  =  r„  +  n,  sec  a  and  n.  is  the 

6  x  0  r  0  o  6 

value  of  boundary-layer  thickness  measured  normal  to  the  body  surface. 

The  values  of  integral  parameters  derived  from  the  measured  axial  velocity 

profiles,  u^U^,  are  given  in  Table  4.  Figure  8  shows  the  theoretical  and 

measured  values  of  (radial)  boundary-layer  thickness  6  and  the  (radial) 

*  r 

displacement  thickness  6  .  The  agreement  between  measured  and  computed 
*  r 

values  of  <5^  and  6^  is  satisfactory. 

MEASURED  TURBULENCE  CHARACTERISTICS 
The  cross-wire  probe  was  used  to  measure  the  turbulence  character¬ 
istics  in  the  thick  boundary  layer.  The  measured  Reynolds  stresses  and 
the  measured  mean  velocity  profiles  were  used  to  obtain  eddy  viscosity  and 
mixing  length. 


MEASURED  REYNOLDS  STRESSES 

The  exchange  of  momentum  in  the  turbulent  mixing  process  is  repre- 

-  2  2 

sented  by  the  distribution  of  the  Reynolds  stresses  -u'v',  u'  ,  v'  ,  and 
2  2  2  2 

w'  where  u'  ,  v'  ,  and  w'  are  the  turbulent  fluctuations  in  the  axial, 
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radial,  and  azimuthal  directions,  respectively.  With  the  probe  elements 
alined  vertically,  the  Reynolds  stress  and  the  axial  and  radial  fluc¬ 
tuations  are  computed  by  an  on-line  computer  with  a  sampling  rating  of 

2 

1024  data  values  in  8  seconds.  To  obtain  the  azimuthal  fluctuation  w'  , 

the  hot  film  sensor  is  rotated  90  degrees.  Data  points  are  accepted  at  the 

same  rate  and  analyzed  by  the  same  computer. 

Figure  9  shows  the  nondimen sionalized  measured  distributions  of 
-  2 

Reynolds  stress  -u'v'/U  and  the  three  components  of  turbulent  fluctuations 


V^7/Un,  >^/Un,  and/w^/ 


V  u'  /U  V v'  /Uq,  andy  w'  /U^  at  various  axial  locations  along  the 
axisymmetric  body.  The  axial  fluctuation  component  is  the  largest  of  the 
three  components  for  all  axial  locations.  The  radial  velocity  fluctuation 
component  has  the  smallest  value.  The  intensity  of  each  of  the  three 
components  is  larger  near  the  surface  of  the  stern  and  reduces  as  the 
boundary-layer  thickens  for  x/L  <  0.95.  A  sharp  gradient  exists  in  the 
turbulence  velocity  components  near  the  surface  since  velocities  go  to 
zero  at  the  wall. 

-  2 

The  maximum  value  of  the  nondimen sionalized  Reynolds  stress  -u'v'/Uq 
occurs  near  the  wall  for  x/L  <  0.873  with  a  strong  variation  in  the  radial 
direction.*  As  the  axial  location  increases  and  the  boundary  layer  thick- 
ends,  the  radial  location  of  maximum  Reynolds  stress  moves  away  from  the 
wall.  The  Reynolds  stress  reduces  quickly  from  the  maximum  value  to  zero 
at  the  edge  of  the  boundary  layer.  At  axial  locations  of  x/L  =  0.951  and 
0.987,  the  measured  values  of  -u'v'  reduce  from  the  maximum  value  to 
zero  at  the  wall.  These  results  are  consistent  with  those  of  Huang  et  al.^ 
for  the  two  afterbody  shapes  tested. 

Huang  et  al.^  also  have  presented  measured  distributions  of  a  turbu- 

-  2  2  2  2  2 

lence  structure  parameter  a^  where  a^  =  -u'v’/q  and  q  =  u'  +  v'  +  w'  . 

Thin  boundary  layer  data  show  a.,  to  have  a  constant  value  of  approximately 

0.15  between  0.056  and  0.86.  The  data  of  Huang  et  al.  indicated  that  for 

a  thick  axisymmetric  boundary  layer,  a^  had  a  value  of  approximately  0.16 

*The  spatial  resolution  of  the  "X"  hot-film  probe  used  may  not  be  fine 
enough  to  measure  the  Reynolds  stresses  precisely  near  the  wall  of  the  thin 
boundary  layer. 
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Measured  Distributions  of  Reynolds  Stresses  for  Afterbody 


Nondimensional  Axial  Lengths 


0.98/ 


up  to  (r-rQ)/6r  =  0.6,  and  that  the  value  of  decreased  toward  the  edge 
of  the  boundary  layer.  The  values  of  measured  on  the  present  axisym- 
metric  body  are  shown  in  Figure  10.  The  present  results  are  similar  to  the 
earlier  results  of  Huang  et  al.^  The  values  of  a^  decrease  in  the  inner 
region  at  the  location  of  x/L  =  0.987. 

No  attempt  was  made  to  remove  the  free-stream  turbulence  fluctuations 
from  the  measured  components  of  turbulence  velocities.  The  reduction  in 

the  value  of  a  near  the  edge  of  the  boundary  layer  is  in  part  caused  by 

1  2  _ 

the  larger  contribution  of  the  free-stream  turbulence  to  q  than  to  -u'v' 

(Reference  1). 


0  0  0.2  0.4  0  6  0.8  1  0  1.2  1  4 


Figure  10  -  Measured  Distributions  of  Turbulent  Structure  Parameter 


EDDY  VISCOSITY  AND  MIXING  LENGTH 


The  experimental  values  of  eddy  viscosity  e  and  mixing  length  £ 

across  the  thick  stern  boundary  layer  are  derived  from  the  measured  values 

_  3u 

of  Reynolds  stress  -u'v'  and  the  mean  velocity  gradient  -5 — ,  using  the 

1  5  13  dr 

definitions  ’’ 


(13) 


ani 

„  3u  9u 

,  ,  „2  x  x  /i#\ 

-u'v'  =  £  37-37-  (14) 

where  the  Reynolds  stress  -u'v'  and  the  axial  velocity  u^  are  measured  by 
the  "X"  hot-film  probe.  By  fairing  the  measured  velocity  profile  ux  with 
a  'spline'  curve,  the  value  of  Su^/3^  can  be  computed  numerically. 

The  experimentally-determined  distributions  of  eddy  viscosity  e/U^  6^ 
are  given  in  Figure  11  for  the  axisymmetric  body  under  consideration.  The 
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Figure  11  -  Measured  Distributions  of  Eddy  Viscosity  for 

Afterbody  5 
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solid  curve  in  this  figure  is  the  thin  boundary  layer  eddy-viscosity 

model  of  Cebeci  and  Smith. The  experimental  values  of  eddy  viscosity  for 

axial  locations,  x/L  0.873,  where  the  boundary  layer  is  thin,  agree 

reasonably  well  with  the  Ccbeci  and  Smith  model.  However,  as  the  boundary 

* 

laver  thickens,  the  measured  values  of  e/U.  6  become  nearly  one-fifth  of 

o  P  5 

the  values  for  thin  boundary  layers  given  by  the  Cebeci  and  Smith  model. 

The  experimentally-determined  distributions  of  mixing  length  P,/5  are 

shown  in  Figure  12.  The  solid  curve  in  this  figure  gives  the  thin  boundary 

13 

layer  results  of  Bradshaw,  Ferrist,  and  Atwell.  Again,  agreement  is 

reasonable  for  x/L  0.873  and  the  experimental  values  in  the  thicker 

boundary  layers  are  roughly  one-third  the  values  in  the  thin  boundary 

layers.  Similar  reductions  of  eddy  viscosity  and  mixing  length  also  were 

1  9 

measured  earlier  by  Huang  et  al.  and  by  Patel  and  Lee.  The 
experimentally-derived  eddy-viscosity  and  mixing-length  data  are  also 
given  in  Table  4. 


Figure  12  -  Measured  Distributions  of  Mixing  Length  for 

Afterbodv  3 


TURBULENCE  MODEL 


Figure  12  shows  that  the  distribution  of  measured  mixing  length  in 

the  thin  boundary  layer  is  roughly  approximated  by  the  thin  boundary  layer 

13 

results  of  Bradshaw  et  al.  However,  the  thick  boundary  layer  values  of 
Jl/6  are  considerably  smaller  than  Bradshaw's  thin  boundary  layer  results. 
Figure  12  also  shows  that  the  maximum  values  of  the  mixing  length  occur  at 
approximately  the  same  radial  location  for  theory  and  experiment.  Huang 
et  al.1  proposed  that  the  mixing  length  of  an  axisymmetric  turbulent 
boundary  layer  is  proportional  to  the  square  root  of  the  entire  turbulence 
annulus  between  the  body  surface  and  the  edge  of  the  boundary  layer 


a 


(VV2"ro 


(15) 


The  present  measured  values  of  (Tq+6  )^-rg  are  given  in  Figure  13 

where  the  solid  curve  is  derived  by  Huang  et  al.1  This  figure  shows  that 
the  values  measured  agree  reasonably  well  with  Huang's  hypothesis. 


Figure  13  -  Similarity  Concept  for  Mixing  Length  of  Thick 
Axisymmetric  Stern  Turbulent  Boundary  Layer 
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r 


The  new  mixing  length  definition  has  been  incorporated  into  the 

3 

Douglas  C-S  boundary  layer  method  of  Wang  and  Huang.  The  original  mixing 
length  formulation  is  approximately  valid  and  is  used  until  the  boundary 
layer  thickness  is  greater  than  20  percent  of  the  body  radius  (6r>0.2  r^) . 
At  this  point,  the  boundary  layer  is  considered  thick  and  the  mixing 
length  i  is  defined  by 


a  =  %, 


/ 


<ro+St)2-ro 


3.31  5 


(16) 


Huang  derived  this  equation  and  the  limits  on  its  use  from  the  maximum 

13 

values  of  mixing  length  given  by  Bradshaw,  £.,  and  the  experimental  data 
1  ° 

of  Huang  et  al. 


TURBULENCE  REYNOLDS  NUMBER 

The  axial  turbulence  velocity  Reynolds  number  is  defined  as 


RA  = 


A 


u'2  A 


(17) 


where  A  is  the  microscale  of  the  turbulence  computed  by  the  relationship 


,2  2 

0  u'  u 

A2  =  - * 

.\2 


(18) 


(?) 


Here,  Taylor's  hypothesis,  the  assumption  of  an  isotropic  field  of  turbu¬ 
lence  (although  not  quite  true  in  Figure  10)  and  space-time  equivalence, 

is  used  to  obtain  the  microscale  A.  The  time  differentiation  of  u'  was 

x 

accomplished  by  using  an  analog  operational  amplifier  (as  used  by  Frenkiel 
14 

et  al.  )  in  a  differentiating  mode.  The  departure  from  linearity  as  a 


fPHPHNPIf 


TTTTJfF 


function  of  frequency  was  found  to  be  1  percent  at  10  kHz  and  5  percent  at 
20  kHz.  The  axial  velocity  fluctuation  u^  was  recorded  on  magnetic  tape 
for  30  sec  at  a  tape  speed  of  152. A  cm/s  using  a  multichannel  recorder. 

This  tape  recorder  had  a  rapid  fall-off  in  response  above  20  kHz.  The 
recorded  data  for  u^  were  played  back  to  determine  the  mean  square  values 
2  - 2 

of  u^  and  (3uV8t)  (with  the  signal  passing  through  the  differentiation 
amplifier,)  by  a  true  RMS  voltmeter. 

The  measured  values  of  turbulence  Reynolds  number  R-^  at  various  axial 
locations  of  Afterbody  5  are  shown  in  Figure  14.  The  value  of  R-^  decreases 


Figure  14  -  Measured  Turbulent  Reynolds  Numbers 


with  increasing  axial  distance  x,  which  is  a  distinguishing  characteristic 

of  a  thick  stern  boundary  layer.  For  a  given  x/L,  the  measured  values  of 

R^  are  almost  constant  for  distances  less  than  60  percent  of  the  boundary 

layer  thickness.  The  measured  values  of  R,  in  this  thick  stern  boundary 

A  14 

layer  are  higher  than  those  generated  by  most  turbulence  grids.  The 
turbulence  Reynolds  number  R^  in  the  present  thick  stern  boundary  layer  is 


i 
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sufficiently  high  that  the  results  obtained  should  be  typical  of  high 
Reynolds  number  thick  stern  turbulent  boundary  layer  flows. 


TURBULENT  LENGTH  SCALE  IN  THICK  STERN  BOUNDARY  LAYER 
The  experimentally  derived  mixing  length  data  in  the  thick  stern 
boundary  layers  of  Afterbodies  1,  2,  and  5  indicate  that  the  mixing-length 
scale  is  proportional  to  the  square  root  of  the  annular  area  between  the 
body  surface  and  the  edge  of  the  boundary  layer.  It  is  of  fundamental 
interest  to  perform  an  independent  and  more  direct  measurement  of  the 
turbulence  length  scale  to  support  this  empirical  result.  The  spatial 
correlation  coefficient  of  two  axial  velocity  fluctuations  at  points  A  and 

B,  p  (d  )  =  u!u'/ (/up^  / u'2) ,  as  a  function  of  the  radial  separation 
AB  r  Ad  A  d 

distance  between  the  points  A  and  B,  =  rg  ”  rA«  is  chosen  for  this 
purpose.  The  axial  velocities  at  points  A  and  B  were  measured  by  position¬ 
ing  two  single-element  TSI  (Model  1212)  hot-wire  probes  with  90-deg  bends 
at  the  two  points.  The  linearized  outputs  of  the  velocities  from  the  two- 
channel  TSI  (Model  1050-1)  hot-wire  anemometer  were  a.c.  coupled  by  using 
an  amplifier  with  low  frequency  rolloff  set  at  1  Hz.  The  fluctuating 
velocity  correlation  u'u'  was  then  measured  by  a  TSI  Correlator  (Model 
1015C).  The  spatial  correlation  coefficients  were  measured  at  various 
axial  locations  in  the  thick  stern  boundary  layer  of  Afterbody  5.  At 

each  axial  location,  point  A  was  fixed  at  r.  =  r_  +  0.167  6  and  the  radial 

AO  r 

location  (r  >r.)  of  point  B  was  increased  until  the  value  of  p._  became 
d  A  Ad 

"small." 

The  measured  spatial  correlation  coefficients  are  plotted  in  Figure  15 


(data  are  tabulated  in  Table  5)  against  three  nondimens ional  length 

/  2  2  * 

parameters,  namely,  n  =  d^//  (r^+d^.)  -r^,  d^/6^.,  anc*  d^/Sp.  As  shown  in 

Figure  15,  the  measured  correlation  coefficients  at  various  axial  locations 

* 

are  not  similar  when  n  is  taken  to  be  d  / 6  or  d  /6  ;  the  correlation 

r  r  r  p 

length  scale  of  the  turbulence  in  the  thick  stern  boundary  layer  is  not 


proportional  to  the  local  radial  boundary-layer  thickness  6  or  the  planar 

*  r 
displacement  thickness  5  .  However,  it  is  interesting  to  note  that  the 
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Figure  15  -  Similarity  Length  Scale  in  Thick  Axisymmetric  Stern 

Boundary  Layer 
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Figure  15b  -  Nondimens ional  ized  by  Planar  Displacement  Thickness 


45 


f 


Figure  15  (Continued) 
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Figure  15c  -  Nondimens ionalized  by  Boundary  Layer  Thickness 


measured  correlation  coefficients  at  various  axial  locations  in  the  thick 
stern  boundary  layer  do  possess  a  similarity  form  when  the  nondimens ional 

/  2  2 

length  parameter  n  is  taken  to  be  d^//  (Tq+S^.)  ~rQ*  This  independent  and 
more  direct  measurement  of  the  turbulence  length  scale  as  determined  from 
the  spatial  correlation  coefficient  measurements  provides  additional  support 
that  the  appropriate  length  scale  in  the  thick  axisymmetric  stern  boundary 
layer  is  proportional  to  the  square  root  of  the  turbulence  annular  area 
between  the  body  surface  and  the  edge  of  the  boundary  layer.  Throughout 
the  entire  stern  boundary  layer,  the  measured  Taylor  microscale  A  at  each 
axial  location  is  found  to  be  approximately  equal  to  the  measured  corre¬ 
lation  length  scale  d^  for  which  the  measured  correlation  coefficient 

p  (d  )  is  equal  to  0.4  (Figure  15).  The  measured  values  of 
AB  r 

A//  (r.+6  )^-r^  were  found  to  vary  from  0.024  to  0.036  for  the  five  axial 
0  r  0 

locations  (Figure  15).  This  is  additional  evidence  of  the  similarity 
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TABLE  5  -  MEASURED  SPATIAL  CORRELATION  FUNCTIONS 


property  of  the  turbulence  in  the  thick  axisymmetric  stern  boundary  layer. 

The  present  set  of  experimental  results  is  consistent  with  the  conjecture 

of  Townsend^  that  the  flow  patterns  in  a  wake  and  in  a  boundary  layer 

possess  an  overall  structural  similarity.  It  may  not  be  necessary  to  use 

"local"  equations  for  turbulent  Reynolds  stresses  with  a  diffusion  term 

9 

and  with  a  correction  to  the  extra  rate  of  strain. 

CONCLUSIONS 

The  present  investigation  of  the  thick  axisymmetric  turbulent 

boundary  layer  on  an  inflected  stern  without  shoulder  separation  is  an 

1  2 

extension  of  the  earlier  research  reported  by  Huang  et  al.  ’  An 

additional  comprehensive  set  of  boundary  layer  measurements,  including 

mean  and  turbulence  velocity  profiles,  turbulence  Reynolds  stresses,  static 

pressure  distribution,  and  correlation  length  scales  of  the  turbulence,  are 

presented.  The  following  major  conclusions  can  be  drawn. 

7  8 

The  Preston  and  Lighthill  displacement  body  concept  has  further  been 

proven  experimentally  to  be  an  efficient  and  accurate  approach  for  computing 

the  viscid  and  inviscid  stern  flow  interaction  on  axisymmetric  bodies. 

The  measured  static  pressure  distributions  on  the  body  and  across  the 

entire  thick  boundary  layer  were  predicted  by  the  displacement  body  method 

to  an  accuracy  within  one  percent  of  dynamic  pressure. 

Neither  the  measured  values  of  eddy  viscosity  nor  the  mixing  length 

were  found  to  be  proportional  to  the  local  displacement  thickness  or  the 

local  boundary-layer  thickness  of  the  thick  axisymmetric  boundary  layer. 

The  measured  mixing  length  of  the  thick  axisymmetric  stern  boundary  layer 

was  again  found  to  be  proportional  to  the  square  root  of  the  area  of  the 

turbulent  annulus  between  the  body  surface  and  the  edge  of  the  boundary 

layer.  This  simple  similarity  hypothesis  for  the  mixing  length  and  the 

displacement  body  concept  have  been  incorporated  into  the  Douglas  C-S 

5  3 

differential  boundary-layer  method  by  Wang  and  Huang.  The  improved 
method  predicts  satisfactorily  the  measured  mean  velocity  distributions  in 
the  entire  stern  boundary  layer. 

The  measured  correlation  length  scales  of  the  turbulence  and  the 
measured  Taylor  microscales  in  the  entire  stern  boundary  layer  are  found  to 


be  proportional  to  the  square  root  of  the  annular  area  between  the  body 
surface  and  the  edge  of  the  turbulent  boundary  layer,  which  indeed  is  the 
appropriate  overall  length  scale  of  the  thick  stern  boundary  layer. 
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1  A.J.  Acosta 
1  I.  Saber sky 
1  D.  Coles 

1  City  College,  Wave  Hill/ 

Pierson 

1  Catholic  U.  of  Amer. /Civil  & 

Mech  Eng 

1  Colorado  State  U./Eng  Res  Cen 


1  U.  of  Connecticut/Scottron 

1  Cornell  U./Shen 

2  Florida  Atlantic  U. 

1  Tech  Lib 
1  S .  Dunne 

2  Harvard  U. 

1  G.  Carrier 

1  Gordon  McKay  Lib 

1  U.  of  Hawaii/Bretschneider 

1  U.  of  Illinois/J.  Robertson 

4  U.  of  Iowa 

1  Lib 

1  L.  Landweber 
1  J .  Kennedy 
1  V.C.  Patel 

3  Johns  Hopkins  U. 

1  0.  Phillips 

1  S.  Corrsin 

1  R.  Long 

1  Kansas  State  U. /Nesmith 

1  U.  of  Kansas/Civil  Eng  Lib 

1  Lehigh  U. /Fritz  Eng  Lab  Lib 

6  MIT 

1  Lib 

1  P.  Mandel 
1  J.R.  Kerwin 
1  P .  Leehey 
1  M.  Abkowitz 
1  J.N.  Newman 

4  U.  of  Minn/St.  Anthony  Falls 

1  Silberman 
1  Lib 
1  Song 
1  R .  Arndt 
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Copies 

4  U.  of  Mich/NAME 

1  Lib 

1  F.  Ogilivie 
1  Hammitt 
1  Cough 

2  U.  of  Notre  Dame 

1  Eng  Lib 
1  Strandhagen 

2  New  York  U. /Courant  Inst 

1  A.  Peters 
1  J.  Stoker 

3  Penn  State 

1  B.R.  Parkin 
1  R.E.  Henderson 
1  ARL  Lib 

1  Princeton  U. /Mellor 

2  U.  of  Rhode  Island 

1  F.M.  White 
1  T.  Kowalski 

1  Science  Application,  Inc. 

McLean,  VA 
C.  Von  Kerczek 

5  SIT 

1  Lib 
1  Breslin 
1  Savitsky 
1  P.W.  Brown 
1  Tsakonas 

1  U.  of  Texas/Arl  Lib 

1  Utah  State  U./Jeppson 

2  Southwest  Res  Inst 

1  Applied  Mech  Rev 
1  Abramson 

3  Stanford  U. 

1  Eng  Lib 

1  R.  Street,  Dept  Civil 
Eng 

1  S.J.  Kline,  Dept 
Mech  Eng 


Copies 

1  Stanford  Res  Inst/Lib 

1  U.  of  Virginia/Aero  Eng  Dept 

1  U.  of  Washington/Arl  Tech 

Lib 

3  VPI 

1  H.L.  Moses,  Dept  Mech 
Eng 

1  D.P.  Telionis,  Dept 
Mech  Eng 

1  J.  Schetz,  Dept  Aero  & 
Ocean  Eng 

3  Webb  Inst 

1  Library 
1  Lewis 
1  Ward 

1  Woods  Hole/Ocean  Eng 

1  Worchester  PI/Tech  Lib 

1  SNAME/Tech  Lib 

1  Bell  Aerospace 

1  Bethlehem  Steel/Sparrows 

Point 

1  Bethlehem  Steel/New  York/Lib 

2  Boeing  Company /Seat tie 

1  Marine  System 
1  P.  Rubbert 

1  Bolt,  Beranek  &  Newman/Lib 

1  Exxon,  NY/Design  Div 

Tank  Dept 

1  Exxon  Math  &  System,  Inc. 

1  General  Dynamics,  EB/ 

Boatwright 

1  Flow  Research 

1  Gibbs  &  Cox/Tech  Info 
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Copies 


Copies 


1 

Grumman  Aerospace  Corp/Lib 

2 

Westinghouse  Electric 

1 

M.S.  Macovsky 

5 

Hydronautics 

1 

Gulino 

1  Lib 

1  E.  Miller 

1  V.  Johnson 

CENTER  DISTRIBUTION 

1  C.C.  Hsu 

1  M.  Tulin 

Copies 

Code 

Name 

1 

Lockheed,  Sunnyvale/Waid 

1 

1500 

W.B.  Morgan 

1 

1504 

V.J.  Monacella 

1 

Lockheed,  California/Lib 

1 

1507 

D.S.  Cieslowski 

1 

Lockheed,  Georgia/Lib 

1 

1508 

R.S.  Rothblum 

McDonnell  Douglas, 

1 

1521 

P .  P ien 

2 

Long  Beach 

1 

1524 

Y.T.  Shen 

1  T.  Cebeci 

1 

1524 

W.C.  Lin 

1  J.  Hess 

1 

1524 

W.  Day 

1 

Newport  News  Shipbuilding/ 
Lib 

1 

1 

1532 

1532 

G.  Dobay 

M.  Wilson 

1 

1541 

P.  Granville 

1 

Nielsen  Eng  &  Res 

1 

1542 

B.  Yim 

1 

Northrop  Corp/Aircraf t  Div. 

1 

1544 

R .  Cumming 

1 

1544 

R.  Boswell 

3 

Rand  Corp 

1 

1544 

E.  Caster 

1  E.R.  Van  Driest 

1 

1544 

S.  Jessup 

1  C.  Gazley 

1 

1544 

K.F.  Lin 

1  J.  Aroesty 

30 

1552 

T . T .  Huang 

Rockwell  International/B . 

1 

1552 

N.C.  Groves 

1 

1 

1552 

G.S.  Belt 

Uj ihara 

1 

1552 

H.T.  Wang 

1 

Sperry  Rand/Tech  Lib 

1 

1 

1552 

1552 

J.  McCarthy 

T.  Langan 

1 

Sun  Shipbuilding/Chief 

1 

1555 

M.S.  Chang 

Naval  Arch 

1 

1560 

G.  Hagen 

1 

Robert  Taggart 

1 

1561 

C . M .  Lee 

TRW  Systems  Group/Lib 

1 

1562 

M.  Martin 

1 

1 

1564 

J.  Feldman 

1 

Tracor 

1 

1568 

G.  Cox 

1 

United  Technology/East 

1 

1572 

M.D.  Ochi 

Hartford,  Conn. 

1 

1572 

E.  Zarnick 

1 

1576 

W.E.  Smith 
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Copies  Code 


Name 


1 

1606 

T.C.  Tai 

1 

1615 

R.J.  Furey 

1 

1802.1 

H.  Lugt 

1 

1802.  2 

F.  Frenkiel 

1 

1840 

J.  Schot 

1 

1843 

H.  Haussling 

1 

19 

M.M.  Sevik 

1 

1940 

J.T.  Shen 

1 

1942 

B.E.  Bowers 

1 

1942 

T.M.  Farabee 

1 

1942 

F.E.  Geib 

1 

1942 

T.C.  Mathews 

1 

1946 

J.A.  Padgett 

10 

5211.1 

Reports  Distribution 

1 

522.1 

Unclassified  Lib 

(C) 

1 

522.2 

Unclassified  Lib 

(A) 
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